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Abstract. A revised (phase 2) single-channel algorithm for aerosol optical thickness, 
retrieval over oceans from radiances in channel 1 (0.63 pm) of the advanced very 
high resolution radiometer (AVHRR) has been implemented at the National Oceanic and 
Atmospheric Administration’s National Environmental Satellite Data and Information 
Service for the NOAA 14 satellite launched December 30, 1994. It is based on careful 
validation of its operational predecessor (phase 1 algorithm), implemented for NOAA 1 1 
in 1989. Both algorithms scale the upward satellite radiances in cloud-free conditions to 
aerosol optical thickness using an updated radiative transfer model of the ocean and 
atmosphere. Application of the phase 2 algorithm to three matchup Sun-photometer and 
satellite data sets, one with NOAA 9 in 1988 and two with NOAA 11 in 1989 and 1991, 
respectively show systematic error is less than 10%, with a random error of <r T — 0.04. 
First results of t£ at retrievals from NOAA 14 using the phase 2 algorithm, and from 
checking its internal consistency, are presented. The potential two-channel (phase 3) 
algorithm for the retrieval of an aerosol size parameter, such as the Junge size distribution 
exponent, by adding either channel 2 (0.83 pm) from the current AVHRR instrument, or 
a 1.6- jam channel to be available on the Tropical Rainfall Measurement Mission and the 
NOAA-KLM satellites by 1997 is under investigation. The possibility of using this 
additional information in the retrieval of a more accurate estimate of aerosol optical 
thickness is being explored. 


1. Introduction 

The upward radiance over dark surfaces, like the open ocean 
in the red and near-infrared spectrum, mainly varies with con- 
centration and type of atmospheric aerosol particles, so that 
information about these particles can be derived from satellite 
data. Aerosols backscatter solar radiation in proportion to the 
aerosol optical thickness i A and the aerosol single scattering 
phase function P A (x ), X being scattering angle. The marine 
aerosol is almost purely scattering (albedo of single scattering 
to -» 1) when not mixed with continental sources (dust, smoke, 
soot) [ Gordon and Morel , 1983; World Climate Program , 1983, 
1986; D 'Almeida et al„ 1991], The relationship between the 
satellite-detected radiance and r* is therefore illumination- 
observation geometry specific. To retrieve A P A (x) must be 
known. 

It was shown in some studies that the variability of P A (x) i n 
backscatter is much less than that of /S Griggs [1975] has 
experimentally shown that upward radiance was essentially 
proportional to t^, independent of aerosol type. Koepke and 
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Quenzel [1979] carried out theoretical calculations for a wide 
range of aerosol models, representative of their natural vari- 
ability, and have shown that variability of P A (x) backscatter 
is well within ^±25% (and can be lowered down to ±4% if 
optimal scattering geometry is used), whereas varies more 
than an order of magnitude. The feasibility of retrieval from 
a single satellite reflectance channel using a “universal” aero- 
sol model has been recently confirmed by Kaufman [1993], who 
has shown experimentally that the relationship between t a and 
atmospheric path radiance is close to universal. Such a one- 
channel procedure may be viewed as a normalization of the 
satellite radiances to a common observation-illumination ge- 
ometry and their scaling to a range of aerosol optical depths, 
measured by customary Sun photometry. The advantage of this 
procedure is that it makes satellite measurements taken in 
different geographical regions and in different seasons uni- 
form, allowing their direct comparison, and scales them to a 
well understandable physical quantity such as the aerosol op- 
tical depth. 

This has been the rationale for the development of an aero- 
sol remote sensing capability at the National Environmental 
Satellite, Data, and Information Sen ice of the National Oce- 
anic and Atmospheric Administration (NOAA/NESDIS), 
which commenced in 1977 [Griggs, 1983]. The algorithm has 
been specifically developed for use with channel 1 (0.63 pm) of 
the advanced very high resolution radiometer (AVHRR) on- 
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AEROSOL OPTICAL THICKNESS(0.6J pm) 

Figure 1 . Dependence of reflected radiation at 0.63 /am on 
aerosol optical thickness and v (Junge particle size parameter) 
as computed from the Dave’ code for the ocean/atmosphere 
model used in the phase 1 retrieval algorithm. Solar zenith 
angle is 60°, satellite zenith angle is 24°, and relative azimuth 
angle is 180°. 


board NOAA polar orbiters, which provides global daily cov- 
erage on an operationally consistent long-term basis. Figure 1 
illustrates the principle of the retrieval of t£ at from the 
AVHRR backscattered radiance. It shows that if the actual 
size of aerosol particles deviates from that used in the retrieval 
model (v ~ 3.5, see section 2 for details), than the 7 ^ AT re- 
trievals would be in error. However, a ground truth experiment 
conducted in 1980-1981 with collocated Voltz-type Sun- 
photometer data sets from 10 island and coastal locations and 
one ship has shown that t£ at retrievals with v = 3.5 have little 
systematic (bias) error, with a random error of ct t ^ 0.05. 

After the phase 1 algorithm was implemented into NOAA 
operations, another ground truth campaign was conducted, 
with emphasis on collocating with higher-quality Sun- 
photometer data from research ships [Reddy et al 1990; 
Kovotyev et al 1993]. The results of this campaign have shown 
a sizeable systematic underestimation by the satellite retrievals 
[Ignatov et al ., 1995a, b]. After taking into consideration the 
poor quality of the 1980-1981 Sun-photometer data, it was 
concluded that the retrieval model used in the phase 1 algo- 
rithm was incorrect and must be revised. The phase 1 algo- 
rithm and results of its validation are described in section 2. In 
its revision, extensive analysis of literature and matchup data 
has been performed (section 3). An improved second- 
generation (phase 2) algorithm has been formulated, validated 
(bias less than 10%; residual random error cr r 0.04), and 
applied to data of NOAA 14, launched December 30, 1994 
(section 4). Both phase 1 and 2 algorithms use only one chan- 
nel, so only one unknown can be retrieved. Even though the 
phase 2 algorithm corresponds to an “average” aerosol better 
than phase 1, it may still be prone to errors in retrieved t£ at 
because of nonuniversality of the aerosol model (resulting 
from seasonal and regional variability in the aerosol physical 
and chemical properties). Measurements in AVHRR channel 
2 (0.83 /am), or, with the launch of NOAA-K in 1996 and the 
Tropical Rainfall Measurement Mission (TRMM) in 1997, a 
third reflectance channel at 1.6 /am, may be potentially used to 
retrieve particle size information, which may be further used to 
reduce errors in t£ at (section 5). The value of the current 
product is discussed in section 6, and concluding remarks are 
given in section 7. 


2. Phase 1 Algorithm (July 1987 to 
September 1994) 

2.1. Algorithm Description 

The phase 1 algorithm for aerosol optical thickness 7 <> AT! 
retrieval has been implemented into an already existing sea 
surface temperature system, which incorporates procedures 
useful for aerosol analysis such as cloud clearing and global 
field analysis [ McClain et al ., 1985]. To avoid the influence of 
specular reflection from the ocean, the solar side of the orbit 
and measurements whose viewing angle is less than 40° away 
from the specular ray have been excluded from the analysis. 
Experimental weekly production of global 8-day (1°) 2 compos- 
ites of 7^ A j j from 2x2 arrays of 4-km global area coverage 
pixels only over oceans commenced in July 1987 [ McClain , 
1989; Rao et al ., 1989; Stowe, 1991]. In July 1989 these global 
analyses began being archived onto magnetic tape at NOAAs 
National Climate Data Center, Asheville, North Carolina. In 
January 1990 they became an operational product at NOAA/ 
NESDIS. 

The principle of the algorithm is shown in Figure I . At the 
particular geometry, the v ~ 3.5 curve is used to retrieve t£ at1 
from cloud-free satellite radiances R at 0.63 pm, reduced to 
dimensionless reflectance by normalization to solar flux, F s , as 
L = tt(R/F s ), and corrected for the Sun-Earth distance. 
T sA , ri » s then scaled to A = 0.50 pm, consistent with the re- 
trieval model (multiplying by ratio of scattering coefficients 
(optical thicknesses) between 0.50 and 0.63 pm, 1.348). In 
operation the relationship between t£ at and L is described by 
a four-dimensional look-up table (LUT), precalculated for dif- 
ferent -E x and view-Sun-azimuth angles using Dave's [1973] 
scalar radiative transfer model, assuming Elterman’s vertical 
profile of aerosol concentration A(h); midlatitude ozone and 
water vapor profiles with integrated contents of 0.316 atm cm 
and 2.91 gm cm 2 , respectively. Other parameters to be spec- 
ified as input to the LUT are Lambertian oceanic reflectance 
p s ; the volume scattering and absorption coefficients, and the 
aerosol phase function, P A {\ ), derived from Mie calculations 
with a prescribed aerosol microphysical model. The oceanic 
albedo is set to p s = 1.5%. The aerosol model uses a modified 
Junge size distribution: dNjdr = 0 (r < r - and r > r V 
dN/dr = A (r min < r < r,„); dN/dr = A(r/r m ) (, ' f {) (r m < 
r ~ ''max)’ where r min , r,„, r max arc particle radii equal to 0.03, 
0.1, and 10 pm, respectively; size parameter v = 3.5 (equiva- 
lent Angstrom exponent a = 1.5); a complex index of refrac- 
tion n — 1. 5—0.0/ (albedo of single scattering, cu, equals 1); 
and A = 22510 cm 3 pm -1 is a normalization constant 
dependent on an assumed number of particles per unit volume 
[Dave, 1973]. This model was adopted from Griggs [1983]. The 
specified modeled LUT establishes a unique relationship be- 
tween the upward radiance and t£ at1 (e.g.. Figure 1), which is 
further used to derive t£ ati globally. 

2.2. Global Application of Phase 1 Algorithm 
to NOAA 11 Data 

The phase 1 algorithm has been applied to data of NOAA 
11, launched in September 1988. Satellite counts have been 
converted to radiances using prelaunch coefficients, with an 
empirical adjustment described by Ignatov et al. [1995a]. This 
product has provided the most complete description to date of 
the spatial and temporal variability of tropospheric and, with 
the eruption of Mount Pinatubo, stratospheric aerosols over 
oceans. Plate 1 shows the variation with time (weekly) from 
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June 29, 1989, to September 8, 1994, of the zonal average (1° 
latitude band) t£ at1 . The most striking feature in Plate 1 is the 
high values in tropical optical thickness coincident with the 
eruption of Mount Pinatubo. The rapid increase and then 
gradual decrease in volcanic aerosol content of the tropical 
stratosphere has been captured, as well as the seasonal trans- 
port of this aerosol to high latitudes [Stowe et al. , 1992; Long 
and Stowe , 1994; Grant et al. , 1996]. Before the eruption and 
beginning about a year after the eruption, other seasonal lat- 
itudinal variations are evident but to a much smaller degree. 
These are associated with natural mechanisms which create 
and disburse aerosols in the troposphere. It is evident from 
Plate 1 that the northern hemisphere has more aerosol than 
the southern. This is due to the aerosol sources primarily being 
restricted to land area (wind-blown dust, fires, and industrial 


pollutants), of which there are more in the northern hemi- 
sphere. 

The geographical separation between tropospheric aerosol 
sources can be seen in Plate 2, which shows the annual mean 
aerosol optical thickness as depicted by this product from the 
2 years before the Mount Pinatubo eruption. It shows that the 
aerosol over the oceans is far from uniform, being separated 
into distinct regional areas of maximum optical thickness. The 
features adjacent to the continents are from wind transported 
aerosol sources originating on the land: mineral aerosol from 
the deserts of Africa, Saudi Arabia, and Asia [Prospero and 
Carlson , 1972]; smoke from biomass burning in central and 
southern Africa, Central and South America, and Indonesia 
[Kaufman et al. , 1990b; Andreae et al. , 1994); and sulfate aero- 
sols from fossil fuel combustion in the highly industrialized 
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Plate 2. Mean distribution of t £ at1 at 0.5 fim for 2 years prior to the eruption of Mount Pinatubo. 


regions of eastern United States, Europe, eastern Asia, and 
Japan [Husar and Wilson , 1993]. There is also a clear indication 
of longitudinal bands of aerosol over the southern and north- 
ern tropical oceans. These are probably the result of biogeni- 
cally produced aerosols over the oceans [Charlson et al. , 1987]. 
Of all the aerosol types present over the oceans, the mineral 
dust is the least similar to the assumed aerosol model and 
therefore is the most uncertain in terms of optical depth esti- 
mates from AVHRR. 

2,3. Validation of the Phase 1 Algorithm 

Sun-photometer data collected during two oceanic 
cruises of the RV Akademik Vernadsky (hereafter AV-89 and 
AV-91, respectively), and one on the NOAA ship Mount 
Mitchell (MM-88; this latter corresponds to NOAA 9 data, 
which were calibrated using prelaunch calibration but without 
the empirical adjustment developed for NOAA 11, referenced 
in section 2.2) were used to quantify the accuracy of the t£ AT] 
retrievals. The principles and results of selecting and process- 
ing the first two matchup data sets are described by Ignatov et 
al. [1995a, b]. The SP measurements taken within ±2 hours of 
the satellite overpasses were selected. The 10 nearest satellite 
retrievals, within 300 km of the ship, were then searched. The 
68 reliable matchups (20 obtained for AV-89, 38 for AV-91, 
and 10 for MM-88) cover large areas in the North Atlantic, in 
different years: 1988 and 1989 (before Mount Pinatubo) and 
1991 (early post Mount Pinatubo), thus being representative of 
a variety of aerosols over the ocean. Averaging of t£ at (over 


the 10 satellite retrievals closest to the ship) and t£ p (averaged 
within ^±2 hours of satellite overpass) was applied to sup- 
press noise resulting from natural variability and measurement 
errors in both data sources. Comparing the averaged charac- 
teristics is also consistent with the gridded NOAA operational 
product which provides an objective analysis of t£ ati fields 
over a ( 1 °) 2 x 1 week space/time box. The accuracy of A r SP is 
-0.02-0.03 [Ignatov et al ., 1995a, b; Reddy et al. , 1990]. 

Figure 2 shows the results of validation of Tsati retrievals 
from NOAA 9 and NOAA 1 1. The regression is 



t A (0.5 ^m) 

SP 

Figure 2. Regression of t£ AX] against Sun-photometer aero- 
sol optical thickness at 0.5-jum wavelength. 
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Plate 3. Global distribution of t £ at2 at 0.63 /am from NOAA 14 (top) for July 1995 with pre-launch 
calibration and (bottom) for August 1995 with postlaunch calibration. 


7s at] = (0.64 ± 0.05) TsV “0.02 (1) 

with only a slightly nonzero bias. The random error of the 
regression is about cr T ~ 0.045, with a correlation coefficient of 


R = 0.87. Thus, to a first approximation, the maps shown in 
section 2.2 should be scaled by the formula (1.56 t£ a , , + 0.03) 
to more correctly represent the optical thickness at 0.50 gm 
observed by Sun photometry. 
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3. Revision of the Phase 1 Algorithm: 

Guidance to Phase 2 

The cause of this bias was the subject of an intense investi- 
gation [Ignatov et al., 1995a]. After removal of the scaling of 
t£ at from 0.63 to 0.50 /am and recalibration of the satellite 
radiances with coefficients determined from postlaunch anal- 
ysis and adjusted for the dates of the matchup data sets [Rao 
and Chen , 1995], the inadequacy of the retrieval model is the 
only reason of discrepancy between t£ at and t£ p . The diag- 
nostic analysis was aided by use of the linearized single scat- 
tering approximation of the radiative transfer equation [Viollier 
et al . , 1980; Gordon and Morel \ 1983], given by the following 
equation for r A retrieval: 
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t sat 4 /Ay /a „ (2) 

where p is an apparent reflectance of the ocean-atmosphere 
system, defined by normalization of L to the solar flux and 
overhead illumination geometry; p R is the Rayleigh scattering 
contribution to the total signal p; p s is the diffuse surface 
reflectance; T is the total atmospheric transmittance; P A and to 
are the aerosol phase function and albedo of single scattering; 

= cos ©„; ja s = cos 0 S ; S u and © s are the view and solar 
zenith angles. Equation (1) shows that the errors in t£ at may 
result from incorrect p s , to, or P A since all other terms are 
reasonably well known. Note that p s participates in (2) as an 
additive term, and (toP A ) as a multiplicative one. This suggests 
that the negative bias in the satellite retrievals (removed in 
operational practice by setting these to zero; note the large 
number of zero values in Figure 2) comes from having over- 
estimated the oceanic reflectance, and the depressed slope 
(^ t sat/^ 7 sp < 1) from an incorrect aerosol model (w and/or 
P A j. 

Analysis of the literature by Ignatov et al. [1995a] suggests 
that p s *** 0.2% rather than 1.5% for deep clean ocean. Cor- 
recting this in the retrieval model and taking into account 
specular reflection of diffuse downward scattered radiation in 
a single scattering approximation [Viollier et al . , 1980; Gordon 
and Morel , 1983] removed the negative bias, but a multiplica- 
tive underestimation in t£ at remained. This result strongly 
indicated that the P A 5 (x) used in the operational retrievals is 
not representative of the average Atlantic aerosol. 

Ignatov et al. [1995a, b] attempted to reconcile t£ at and r£ P 
by inferring a microphysical model from the wavelength de- 
pendence of the t£ p. This gave an Angstrom exponent of about 
a — 0.6, which led to setting v 2.5 (v = a + 2) in the Junge 
size distribution. However, using this model removed the sys- 
tematic discrepancy between Ts AT and t£ p only when it was 
supplemented by adjusting the imaginary part of the aerosol 
refractive index to 0.01. Correlation statistics indicate that this 
algorithm estimates t£ ax with little systematic bias or offset: 
t£ at -0.004 + (0.98 ± 0.08)t£ p . The correlation coefficient 
for the 68 cases is R = 0.92, with a random error of a r ** 
0.046. Additional analyses of Mie calculations have shown that 
the imaginary part influences the results of t£ at retrieval 
mainly through decreasing the phase function in the backscat- 
ter direction (about 30%) rather than through absorption itself 
(only 10%, consistent with the corresponding albedo of single 
scattering, a> *** 0.9). 

The value of a) ~ 0.9 was thought of as unrealistic, and A. 
Ignatov (submitted manuscript, 1996) tried to estimate an em- 
pirical phase function, P A m pJ (Xj), at each matchup point i for 


Figure 3. Two model phase functions: P A S (Junge, v = 3.5) 
(dashed) and P A N derived from Mie calculations with lognor- 
mal size distribution (r m = 0.1 pm; o = 2.03 jam; n = 1.40- 
0.0/] (dotted), and results of the empirical phase function, 
F A mp , estimation for AV-89 (crosses) and AV-91 (open cir- 
cles). All phase functions are normalized to 477. 


the particular scattering angle Xv having the r^ P known from 
the Sun-photometer measurements. The retrieved values of 
PempAxd are plotted versus scattering angle in Figure 3, each 
point corresponding to one matchup. Ail the experimental 
points fall well below P A 5 , which explains the systematic un- 
derestimation in the operational retrievals (Figure 2 and (2)). 
One has therefore to find an aerosol phase function which 
describes the experimental points in Figure 3. 

According to D’Almeida et al. [1991], clean maritime aerosol 
is composed mostly of two components: sea salt (in which the 
nucleation mode dominates) and sulfate. For a typical relative 
humidity over the oceans of 70-80% the lognormal distribu- 
tions of each are close together with r m «■ 0.1 /am, and cr - 2.03 
jam. Complex indices of refraction for both are also close, 
about 1. 4-0.0/. As a result, the bimodal size distribution of 
the marine aerosol can be effectively represented by a mono- 
modal lognormal function. Its phase function, shown in Figure 
3, describes the experimental points quite well. Figure 4 shows 
7 ^ at retrieved with the described lognormal model versus 
A r SP . The regression gives 

t&t - 0.01 + (0.91 ± 0.04)4. (3) 

with R — 0.93 and an error of cr T 0.04. This error is inherent 
to the present single-channel algorithm, and can possibly be 
reduced with future multispectral algorithms. 

Note that the A. Ignatov (submitted manuscript, 1996) study 
includes only matchups with large aerosol loadings (t£ p >: 0.1), 
which are not typical for clean maritime conditions [Smirnov et 
al ., 1995]. As a result, the derived empirical phase function 
may not be representative of clean maritime aerosol, and one 
needs to be cautious in the geophysical and microphysical 
interpretation of the derived phase function. Also, application 
of the derived phase function in the case of small aerosol 
loadings may be questionable. However, for practical retrievals 
from AVHRR, error in P A is of more concern for large aerosol 
loadings, since absolute systematic error <5 t£ at is proportional 
to Tsat f° r the same relative error in P A (see equation (2)). 
Note that reducing the absolute uncertainty 8 t£ at is important 
for climate modeling [e.g., Charlson et al ., 1992]. To that end, 
the lognormal phase function is a reasonable candidate for 
replacing, the P A 5 (x) used in the phase 1 retrieval algorithm. 
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Figure 4. Regression of t£ AT 2 against same Sun-photometer 
data set as in Figure 4 but at 0.63 pm. 


The phase 2 model's apparent contradiction to other litera- 
ture, which recommends using a two- or three-mode lognormal 
function [ World Climate Program, 1983, 1986; Shettle and benn. 
1979; Gathman, 1983; Gras, 1991; Hoppelet al.. 1994], requires 
further investigation hut is of little concern for practical use, 
since the P, X N is close to the empirically derived one. The 
premise of that model is its simplicity, in particular fewer 
numbers of descriptive parameters (each mode requires values 
of modal radius, standard deviation, and complex index of 
refraction). This model has been adopted tor the operation ol 
NOAA 14, launched in December 1994. Its summary and some 
preliminary attempts at validation are described in the next 
section. 

4. Phase 2 Algorithm (March 1995 Onward) 

4.1. Algorithm Summary 

The second-generation (phase 2) aerosol retrieval algorithm 
has many features common with phase 1 described in section 2. 
The difference between the two comes from the fact that the 
phase 1 algorithm has been carefully validated, and its param- 
eters have been revised based on experimental data, and care- 
ful analysis of the recent literature. The Lambertian surface 
reflectance in the second generation product is lowered 7.5 
times, from 1.5% down to 0.2%. The bidirectional Fresnel 
reflectance from the flat ocean surface is approximately taken 
into account by introducing a diffuse glint correction to the 
phase function, derived from the single scattering approxima- 
tion [Vi oilier et al 1980; Gordon and MoreL 1983]. Its appli- 
cation to the phase 2 algorithm is described in detail by Ignatov 
et al. [1995a]. We understand restrictions of this procedure, 
especially at high aerosol loadings and in the vicinity of glint, 
and are working now on more appropriate quantification of 
this effect. The aerosol phase function used in the phase 2 
algorithm, was derived assuming a lognormal size distribution 
with r m = 0.1 /am, a = 2.03 /am; n = 1.4-0/. 

No scaling is made anymore from the wavelength of mea- 
surement 0.63 /am to the wavelength of 0.50 /am, recom- 
mended by the World Meteorological Organization as a refer- 
ence. The reason is that this scaling factor is strongly 
atmospheric model dependent and may result in a systematic 
error in the retrieved t£ a , . For those wanting to convert, the 
scaling factor for the first-generation model is 1 .348 and for the 
second-generation is 1.164. 


4.2. Testing Phase 2 Algorithm With NOAA 11 Pathfinder 
Data (September 1989) 

Expected differences in the performance of the two algo- 
rithms has been verified globally with data from the AVFIRR 
Atmosphere Pathfinder Program. This is a NASA/NOAA joint 
program to reprocess NOAA satellite data with community 
consensus calibration and remote sensing algorithms into a 
long-term product data set for weather and climate change 
research [University Corporation for Atmospheric Research 
( UCAR ), 1994]. The AVHRR atmosphere pathfinder data set 
is being produced at NOAA/NESDIS. Phase A l of this pro- 
gram uses the CLAVR-1 (clouds from AVHRR phase 1) cloud 
screening algorithm [Stowe et al .. 1991], to produce among 
some 71 parameters, cloud-free and cloudy radiances in all five 
channels averaged in 1 10 km equal area regions over the globe. 
Phase A2 uses these area averaged radiances to derive aerosol 
optical thickness over oceans and radiation budget parameters 
at the top of the atmosphere. Initially, an 18-month “bench- 
mark’' period from April 1987 to November 1988 (NOAA 9) 
and September 1989 (NOAA 11) is being processed. For this 
initial run, t£ at has been computed with both phase 1 and 2 
aerosol algorithms for validation purposes and for continuity 
with the over 5 years of NOAA 1 1 data. Figure 5 show's the 
result of plotting t£ at1 versus t£ AT 2 (converted to 0.50 pm for 
consistency) for September 6, 1989, NOAA 11 data tor all 
1 10-km grid cells with cloud-free radiances. The regression line 
has the form 

Tsai , = (0.69 ± 0.003 )t£ VI2 - 0.03 (4) 

with R - 0.97 and <r T - 0.023. The relationship between t£ ati 
and t£ p (cf. equation (1 )) produced a slope and intercept which 
agrees to those in (3) within the error of estimation. Thus it can 
be concluded that when globally applied, the phase 2 algorithm 
yields a result that is similar to what would have been observed 
had we performed validation with Sun photometers. 

4.3. Application of the Phase 2 Algorithm to NOAA 14 
Data 

After the failure of NOAA 1 1 AVHRR in September 1994, 
operational production of t£ AT 2 did not commence until 
March 1995 using NOAA 14 data. As an example of these 
initial results, Plate 3 shows a global image of t£ AI2 for the 
month of July 1995. It reproduces many of the features present 



t A (0.5 fim) 

SAT2 

Figure 5. Regression of t£ at , against t£ AT2 at 0.5 pm glo- 
bally, using the AVHRR atmosphere pathfinder processing 
system for September 6, 1989. 
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Figure 6. Effect of prelaunch calibration error on t£ AT2 . 


in the NOAA 1 1 maps (cf. Plate 2), but the optical depths 
appear to be elevated by more than the expected factor of 1 .34 
(a ratio of underestimation of 1.56 from (1) and conversion 
factor from 0.63 to 0.5 /am of 1.164). Postlaunch investigations 
into the calibration of NOAA 14 AVHRR channels 1 and 2 (N. 
Rao, personal communication, 1995), have confirmed revision 
of prelaunch calibration by R. Mitchell (In-flight characteris- 
tics of NOAA AVHRR channels 1 and 2, submitted to IEEE 
Transactions on Geoscience and Remote Sensing , 1996) and 
have led to the conclusion that the prelaunch values of cali- 
bration coefficients were in error. Rao’s revision of these terms 
was implemented into operations the end of July. The effect on 
the NOAA 14 aerosol product is shown also in Plate 3, where 
much more reasonable values of t£ AT 2 appear to be distributed 
over the globe for August 1995, having about the right degree 
of enhancement compared with Plate 2. Figure 6 shows quan- 
titatively the effect of applying Rao’s postlaunch calibration 
coefficients to a sample of NOAA 14 data. This shows that 
t£ AT2 was being overestimated before August 1995 by about 
0.08, on average. 

The Sun-photometer method of validation is nearing a re- 
ality for NOAA 14. NESDIS has begun collecting Ts AT 2 ob- 
servations in real time at seven coastal locations where well 
calibrated and maintained Sun-photometer observations are 
available, but none of these data sets have been processed as 
yet [Holben et ai, 1996]. 

5. Possible Enhancements to the AVHRR 
Product 

5.1. Multispectral Retrieval Algorithms 

A third-generation (phase 3) algorithm is currently under 
development. It uses two channels from the current AVHRR. 
The algorithm will be used to estimate an additional parame- 
ter: the Angstrom wavelength exponent a traditionally mea- 
sured in Sun photometry. This could be used to retrieve an 
aerosol size parameter, which in turn could possibly be used to 
adjust the size distribution (scattering phase function) used in 
the single-channel retrieval algorithm [e.g., Durkee et ai, 1991]. 
This could potentially remove residual errors resulting from 
departures of the actual aerosol size distribution from the one 
used in the single-channel retrieval model. 

One example of such an algorithm was developed several 
years ago [Ahmad et al., 1989]. This approach relied upon 
regression relationships from theoretical radiative transfer 


model calculations which relate the ratio of aerosol path radi- 
ance between two channels to the exponent of the Junge power 
law model. The retrieved Junge exponent (v = a + 2) deter- 
mines the aerosol particle size distribution. It is then used to 
retrieve the optical thickness from one of the channels, again 
through the use of theoretical regression relationships, but 
specifically for that particular Junge exponent (aerosol size 
distribution). 

A graphical illustration of such an algorithm is shown in 
Figure 7a. On the ordinate normalized radiance (defined in 
section 2.1) is plotted at the effective wavelength (average 
wavelength weighted by product of channel response function 
and solar spectral irradiance) of channel 2 (0.83 /am) of 
AVHRR, and on the abscissa the same quantity is plotted but 
for channel 1. Solar and satellite geometries are identical to 
those used in Figure 1. This has been done for three atmo- 
spheric models (tropical, midlatitude, and Arctic) to represent 
the possible range of water vapor content which affects the 
radiance in channel 2 due to its filter encompassing several 
strong bands of water vapor absorption. To each of these three 
atmospheric models is added the ocean surface reflection 
model and Junge aerosol size distribution dependent scatter- 
ing and absorption models used in our initial efforts to improve 
upon the phase 1 single-channel algorithm [Ignatov et al., 
1995a]. Three values of the Junge exponent have been used in 
this example: v = 3.5 (effective radius, third moment divided 
by second moment of size distribution, r cfi = 0.21 /urn), v = 2.5 
(r cff = 0.79 jum) and v = 1.5 (r cff = 3.64 /urn). The total 
number of particles in each model atmosphere was adjusted to 
yield the same aerosol optical thickness at 0.5 /urn (heavy 
dashed line in Figure 7 shown only for midlatitude model). The 
result is a set of three families of curves, one for each Junge 
exponent value, relating the change in radiance for the two 
channels as optical thickness changes. Making up each of these 
are three curves for each of the three different atmospheric 
models, showing the effect of variable water vapor on the 
radiance in channel 2. 

Given an AVHRR observation over a clear ocean at the 
viewing and illumination angles of Figure 7 (it should fall 
somewhere in the space bounded by the outermost curves), the 
corresponding values of v and t£ AT 3 can then be read. Auto- 
mating this procedure is relatively trivial, once a similar set of 
figures (i.e., look-up tables) have been computed for the range 
of geometries likely to be observed. The accuracy of the re- 
trieved optical thickness and Junge parameter is dependent 
not only on the correctness of the assumed model parameters 
but also on measurement error. As is evident from Figure 10a, 
measurement noise will introduce larger relative error in the 
derived parameters as the optical thickness of the aerosol ap- 
proaches zero (curves converge at t£ AX 3 = 0, where the ability 
to separate the effect of different aerosol models vanishes). 
Quantifying the magnitude of these errors and their sensitivity 
to independent variables is the subject of current research 
(e.g., A. Ignatov, Accuracy of aerosol Angstrom exponent de- 
rived from two-channel satellite measurement over the oceans, 
submitted to Applied Optics , 1996). 

One source of error, the variable amount of water vapor, can 
be reduced in two ways: either by providing independent 
knowledge of the water vapor amount and vertical profile 
relative to the aerosol profile, not presently available from 
satellite remote sensing technology, or by using a pair of chan- 
nels weakly affected by water vapor absorption. Such a pair of 
channels (0.63 and 1.61 /urn) will become available from 
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Figure 7. Illustration of the principle of t£ AT 3 and v retrieval from two-channel reflectance measurements 
at (a) 0.63 and 0.83 pirn and (b) 0.63 and L61 jam (the heavy dashed lines are contours of t£ at3 at 0.5 pjn, 
shown only for the midlatitude model). 


AVHRR with the NOAA-KLM series of satellites, scheduled 
to commence in 1996. However, because of data storage con- 
straints on-board the spacecraft, only the afternoon satellite, 
NOAA-L, will provide global coverage at 1.61 /xm in daytime. 
To show the potential of retrieving improved aerosol informa- 
tion from these two channels. Figure 7b has been constructed. 
Not only is the effect of variable water vapor greatly reduced, 


but the separation of the three families of curves has been 
increased. This means that for the same relative measurement 
error, one should be able to obtain more accurate aerosol 
information at lower values of optical thickness than with the 
current two AVHRR reflectance channels. A similar pair of 
channels are available from the Thematic Mapper on the 
Landsat series of satellites, from the along track scanning ra- 
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diometer (ATSR) instrument on the ERS-2 satellite, and from 
the visible infrared scanner (VI RS) instrument to be launched 
on the TRMM satellite in 1997. 

5.2. Other Potential Enhancements 

With regard to other aerosol remote sensing issues not yet 
addressed in our program, it is highly desirable to make re- 
trievals over surfaces with relatively high reflectivity, such as in 
the regions of oceanic specular reflection or for most land 
surfaces. A climatology of clear-sky reflectance and its variance 
could be used as input to radiative transfer model LUTs for the 
extraction of optical thickness and particle size parameters in 
these areas, although the sensitivity of reflected radiation to 
the presence of aerosols is diminished relative to that over dark 
ocean surfaces [ Kaufman et al., 1990a]. Also, use of infrared 
retrieval techniques, as other investigators have shown [Leg- 
rand et al ., 1989], can retrieve aerosol optical thickness infor- 
mation over certain geographical areas, such as deserts. Even- 
tually, with the launch of the ocean color temperature scanner 
(OCTS) and polarization and directionality of Earth reflec- 
tances (POLDER) instruments on the Japanese ADEOS sat- 
ellite and SeaWifs instrument on a commercial satellite in 
1996, and the launch of the EOS polar platform instruments 
(moderate-resolution imaging spectrometer (MODIS), mul- 
tiangle imaging spectrometer (MISR), and Earth observing 
scanning polarimeter (EOSP)) between 1998 and 2003, many 
other spectral channels as well as linear polarization parame- 
ters will become available for the description of the atmo- 
spheric aerosol content. 

Also, even at the single-channel algorithm level, further im- 
provements could be made by more accurately describing the 
diffuse glint and introducing the direct glint, both of which will 
require information about wind speed. Using a microphysical 
aerosol model that varies by geographical region and provides 
a phase function that closely matches empirically derived val- 
ues for each region separately may also improve accuracy. 
Establishment of ground truth sites with in situ as well as 
remote sensing data would greatly assist the validation of mul- 
lispeetral techniques. 

6. Discussion 

The phase 2 retrievals show reasonable agreement with Sun- 
photometer data on the average. Residual scatter indicates 
that the phase function for real maritime aerosol varies and 
can be assumed “universal” only within the reported uncer- 
tainty limits. For single-channel retrievals the only alternative 
to a “universal” empirical phase function would involve its 
parameterization versus t*: P a = /(r A ), which allows one to 
keep one unknown under retrieval The analysis of A. Ig- 
natov (submitted manuscript, 1996) did not confirm this hy- 
pothesis, however. Also, oceanic diffuse reflectance is not a 
universal constant, p s = 0.2%. It may be much higher over 
turbid coastal waters, and is bidirectional. 

We believe that a more complicated tw'o-channel procedure 
should be introduced only after the potential of the one- 
channel algorithm is fully realized, its advantages and deficien- 
cies understood, and the necessity and sufficiency of additional 
measurements convincingly proved. To that end, using the 
microphysical model which closely reproduces the empirically 
derived phase function is the logical next step, providing a 
bridgehead between the present and future generation of aero- 
sol retrieval algorithms. Additionally, the spectral position of 


AVHRR channel 1 is not far from the optimal for aerosol 
retrieval, recommended by Koepke and Quenzel [1981]. So far, 
the rms accuracy of the Phase 2 algorithm is rr T - 0.04. This 
estimate is consistent with results of Kaufman [1993], who 
shows that using a universal relationship between aerosol op- 
tical depth and path radiance leads to a maximum error of 0.10 
(about 3 (j t , assuming a Gaussian distribution of error). 

It is recognized that the AVHRR instrument is not the 
optimum instrument for making aerosol measurements. It has 
no onboard calibration for its two (a third will be added in 
1997) reflectance measuring channels. With the development 
of spaceborne multispectral imaging spectrometers like MO- 
DIS [King et al ., 1992], multiangular viewing radiometers like 
MISR [Martonchik and Diner , 1992], and ATSR [Mutlow et al. , 
1994], and polarimetric imaging radiometers like POLDER 
[Deitze et al. , 1993], remote sensing of aerosol physical and 
chemical properties should be improved beyond AVHRR’s 
capabilities. Nonetheless, the work that has been done with 
AVHRR has clearly demonstrated that aerosol remote sensing 
can be done from space, and that the resulting products over 
global oceans can provide valuable information on spatial and 
temporal variability of this important climate forcing parame- 
ter [J ous sail me and Sadourny, 1989; Long and Stowe , 1994; 
Lacis and Mishechenko , 1995]. Efforts are now underway to 
reprocess the entire AVHRR record from the NOAA after- 
noon satellites back to 1981 [ UCAR , 1994]. Such a climatology, 
having accounted for calibration drift and the same aerosol 
retrieval algorithm throughout, wall provide a unique data set 
from which a greater knowledge of the Earth’s aerosol char- 
acteristics will be gleaned, and against which the value of the 
future instruments can be measured. 

7. Conclusion 

The second generation single-channel aerosol retrieval al- 
gorithm has been described. It has been implemented into the 
operational processing system for the NOAA 14 satellite, 
launched on December 30, 1994. The algorithm is based on 
analysis of three matchup Sun-photometer-NOAA satellite 
data sets, which revealed errors in the operational NOAA/ 
NESDIS t£ ai1 retrievals. A fundamental revision of the oce- 
anic reflectance model was undertaken from analysis of the 
scientific literature. A new empirical aerosol phase function 
has been derived from the matchup data. A microphysical 
aerosol model, which matches the empirical phase function, 
has been chosen for the retrieval of t£ at . These adjustments to 
the first generation NOAA/NESDIS aerosol product bring, to 
a first approximation, the Sun-photometer and satellite data 
into better than 10% agreement. It also has been shown that 
when applied to a global analysis with AVHRR atmosphere 
pathfinder data sets, the phase 2 algorithm provides r£ An 
which deviates from t£ a1 , in the same way that deviated 
from t£ ai j. Examples of the NOAA 14 aerosol global maps 
have been presented, which before August 1995, are biased 
high by AVHRR calibration errors. 

Since the t£ at are retrieval model dependent, we call the 
satellite product a “radiatively equivalent” aerosol optical 
thickness. The phase 2 algorithm provides unbiased, in the 
mean statistical sense, aerosol retrievals over oceans. It can be 
viewed as a procedure for removing the strong angular depen- 
dencies from upward radiances (normalizing them to a com- 
mon geometry), with scaling to aerosol optical thickness mea- 
sured by Sun photometry. Any remaining noise in the satellite 
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retrievals (cr T » 0.04, see (3)) follows from variability of the 
aerosol phase function, aerosol absorption, and oceanic reflec- 
tance about their average values, typical for the oceans. Some 
of this variability may be potentially accounted for by further 
adding other spectral or angular measurements to the retrieval 
algorithms; for those, the present NOAA/NESDIS aerosol 
product provides a reference level from which improvements 
can be measured. An example of such an algorithm, with the 
improvements possible with the 1.61-p,m channel on the 
NOAA/KLM series of AVHRR instruments, was presented. 
Further investigations are needed, together with intensive field 
measurements, to establish the most globally robust multichan- 
nel aerosol retrieval algorithm. 
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